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Abstract 



\ We prove that any spectral sequence obeying a certain growth 

^sO ' law is the quantum spectrum of an equivalence class of classically in- 

tegrable non- linear oscillators. This implies that exceptions to the 
\ Berry- Tabor rule for the distribution of quantum energy gaps of clas- 

sically integrable systems, are far more numerous than previously be- 
£>v lieved. In particular we show that for each finite dimension k, there are 

an infinite number of classically integrable A;-dimensional non-linear 
oscillators whose quantum spectrum reproduces the imaginary part of 
zeros on the critical line of the Riemann zeta function. 

o 



PCACS: Numbers 05.45, 03.65, 02.30.Dk, 02.30.Gp, 02.10.Lh 



c3 1 Introduction 



An important theme in the quantum theory of classically chaotic systems 
is the relationship between the qualitative behaviour of the classical system 
and statistical properties of its quantum mechanical spectrum A much 
studied statistic in this regard is the distribution of energy level spacings. 
According to a result of Berry and Tabor the values of the energy gaps of a 
generic integrable system are Poisson distributed 0. In general this differs 
remarkably from the statistical distribution of energy gaps in the case of 
classically chaotic systems. These have been studied numerically and can be 
described by Wigner, GOE or GUE rather than Poisson statistics ffl]. 
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An exception to the result of Berry and Tabor occurs in the case of the 
harmonic oscillator. However Razavy investigated a family of integrable per- 
turbations of the harmonic oscillator and found that this departure from 
Poisson statistics is non-generic in the sense that perturbing away from the 
harmonic oscillator, the energy gaps quickly become Poisson distributed |J. 
Another exception to the Berry-Tabor result was pointed out by Casati, 
Chirikov and Guarneri, and occurs for the case of a free particle in a rect- 
angular well 0| . Nevertheless Seligman and Verbaarshot observe that for 
potential wells close to the rectangular well, the distribution of energy gaps 
is close to Poisson. They conclude that in this case also departure from 
Poisson statistics is non-generic ||. The general perception is that despite 
exceptional cases such as the harmonic oscillator and the free particle in a 
rectangular well, the statistics of quantum energy gaps of classically inte- 
grable systems are universally described by the Poisson distribution. 

In this paper we demonstrate that given a spectral sequence obeying a 
certain growth law, there exists an infinite family of classically integrable 
Hamiltonians whose quantum spectrum coincides with this sequence. This 
shows that a wide range of exceptions to the Berry- Tabor rule are possible 
and that any quantum system whose spectrum obeys such a growth law, can 
be simulated by a family of classically integrable non-linear oscillators. In 
particular we derive a result concerning the hypothesis of Berry about the 
qualitative behaviour of an unknown classical system whose quantum energy 
levels are given by the imaginary part of the zeros on the critical line of the 
Riemann-zeta function ||. Berry argues that this classical system should 
be chaotic. We prove that this unknown classical system is not unique and 
need not be chaotic. In fact we show that an infinite number of classically 
integrable non-linear oscillators are capable of reproducing these zeros when 
quantised. 

2 One dimensional case 

To justify our claims we require a number of theorems based on the following 
lemma and its generalisation. 

LEMMA 1 Given a sequence of complex numbers {£„ : < n G Z} obeying 
a growth law |£ n | < exp(a + bn) for some a G 3ft, b G 9ft + , there exists an 
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equivalence class S of entire functions of the complex plane such that for each 
s G S we have s(n) = ( n . 



We prove this by considering 



oo 



(nfn(z) 



(1) 



n=0 



where for 0<nG2, < e G and z 7^ n, 



f n (z) = exp((z -n)(2n + b + e)) 



sin(27r(,2 — n)) 



2tt(z — n) 



Strictly speaking f n has a singularity at z = n. The singularity is however 
removable and so defining f n (n) — 1, f n becomes an entire function of the 
complex plane. Eventually we will show that s is well defined by the series in 
eqn(|l]), and that it too is an entire function of the complex plane. Observing 
that f n (m) = 5 nm for < n, m G Z, the property s(n) = ( n follows from 
the definition of f n and an explicit evaluation of the sum at z — n. This 
construction provides a single repesentative member of the equivalence class 
S. The difference between any two representatives is an entire function of 
the complex plane which vanishes on the non-negative integers. This set is 
infinite and denoting it by 5*0, the equivalence class is given by S = s + S . 
All that now remains is to show that s is entire. 

The domain {z G C : \z\ < p G 3?} will be denoted D p . Using A(D p ) 
to represent the analytic functions on D p , the entire functions are naturally 
denoted A(C). The Weierstrass M-test provides a criterion for when the 
infinite sum of functions which are analytic on some domain D, converges 
to a function which is analytic on D. Specifically if {g n } is a sequence in 
A(Dp), and if there exists a sequence of positive real numbers {A n } with 
the property A n > \\g n \\ = sup^ eDp |g n (z)|, such that Y^=qA u < 00, then 
Y^=o9n converges to g G A(D p ). To prove that s is analytic on an arbitrary 
D p , we apply this test to the sequence defining s in eqn(|l]). Each term lies 
in A(D p ) since each is a constant multiple of an entire function. Using the 
fact that \sin(z)/z\ < exp(|z|) we have for z G D p 



\Un(z)\ < |Cn| |exp((z-n)(27r + & + e))| 



sin(27r(,2 — n)) 



2it(z — n) 
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< exp(a + 6n) exp((|z| — n)(2ir + b + e)) exp(2ir\z — n\) 

< exp(a + bn) exp((|z| — n)(b + 2ir + e)) exp(27i(\z\ + n)) 
= exp(a + 1^1(4^ + b + e)) exp(— en) 

< exp(a + p(4iT + b + e)) exp(— en). 

To apply the Weierstrass M-test we take A n = c p exp(— en) where the con- 
stant c p = exp(a + p(4ir + b + e). The sum X^oAi converges and so 
s G A(D P ). Since s G A(D p ) for any p G 3? + we have s G A{C) and 
lemma [I] is proven. 

Although there is a different s for each valid choice of a, b and e, the 
corresponding equivalence class S is independent of these values. Lemma [I] 
leads almost directly to the following theorem. 

THEOREM 1 Given a real sequence {E n : < n G Z} which obeys a 
growth law \E n \ < exp(a + bn) for some a, b G 3? + , there exists an equivalence 
class of classically integrable non- linear oscillators He such that if h is the 
Hamiltonian of the simple harmonic oscillator, each he G He is of the form 
hc(h°) and its quantum spectral sequence is given by E n . 

To see this construct s as in lemma |T] so that s(n) = E n for < n G Z. 
The equivalence class of classical Hamiltonians He comes from replacing z 
with h° = (q 2 +p 2 )/2 in each representative s G S. As h° is the Hamiltonian 
of the simple harmonic oscillator, each he G He is the Hamiltonian of a 
classically integrable non-linear oscillator. 

Since s G A(C), each he G He can be identified with its Taylor ex- 
pansion in h which is everywhere convergent. This allows us to write 
he = Do° c i (h ) 1 where each q G 3?, and to define a corresponding quantum 
Hamiltonian operator 

CO 

Hq = £ c ^ (2) 
o 

where iV* is the product of z-copies of the number operator. With h~ — 1 
is essentially the Hamiltonian of the quantum harmonic oscillator. This 
procedure quantises the classical Hamiltonian he G He and determines an 
appropriate operator ordering. The action of Hq on L 2 (x, dx) is well defined, 
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the eigenfunctions of Hq are the familiar harmonic oscillator eigenfunctions, 
and the corresponding eigenvalues are simply E n . 

Although the coefficients of Hq in (|2]) depend on the choice of he G He, 
only values of he{z) when z is a positive integer play a role in the dynamics. 
It does not matter which h c G H c we use to construct Hq, the resulting 
quantum system will be the same. Hq is a unique quantum Hamiltonian, 
uniquely determined by the spectral sequence, and which we can identify 
with the equivalence class of classical Hamiltonians He- 

This is a little surprising since it indicates that the correspondence of 
classical to quantum systems is not one-to-many as we might naievely have 
expected. The usual correspondence between classical observables Oc{q,p), 
and quantum observables Oq{q,p,K) is one-to-many in the sense that for 
a given classical Hamiltonian he{q,p), there is an infinite number of cor- 
responding quantum Hamiltonians hg(q,p,h) determined by the property 
hg(q,p,0) = he(q,p)- In our construction the redundancy due to different 
choices of operator ordering does not arise. A different redundancy however 
does arise. This is because there are different classical systems he which 
correspond to the single quantum system Hq. This is directly attributable 
to the fact that many different continuous functions interpolate a function 
whose values are specified only on the integers. For a given value of Plane's 
constant the differences between these classical systems occur on a scale 
smaller than h. For compact dynamical systems such as those used in our 
construction, the correspondence between classical and quantum systems is 
in fact many-to-many. 

3 Finite dimensional case 

It is straightforward to generalise the results of the previous section to the 
k- dimensional case. A natural generalisation of lemma [I] is given by 

LEMMA 2 Given a k-indexed sequence of complex numbers {Cm...n fc : < 
Tii G Z, 1 < i < k} obeying a growth law |Cm...nJ < exp(a + b\ni + . . . + bkUk) 
for a G 3?, h{ G 9ft + , there exists an equivalence class S of entire analytic 
functions on C k such that s{ri\, . . . , nk) = (ni...n fe for all s G S. 

This leads to a natural generalisation of theorem [I] as follows. 
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THEOREM 2 Given a k-indexed sequence of real numbers {E ni ... nk : < 
rii G Z , 1 < i < k}, which obeys the growth law \E ni ___ nk \ < exp(a+b\ni + . . .+ 
bk^k) for some a G 3ft, bi G 3ft + , there exists an equivalence class of classically 
integrable k-dimensional non-linear oscillators He, where if h® for 1 < % < k 
are independent classical harmonic oscillator Hamiltonians, each h c G He 
is of the form h c (hi, . . . , hi) and its quantum spectral sequence is i? ni ... nfe . 

The proofs of lemma |2] and theorem |2| are based on a consideration of 

oo 

s(z\, . . . , Zk) — ^ (ni...n k /ni(^l) - • • • " fn k ( z k) 
rai...ra fc =0 

where for < q G 9ft we have 

r- i \ a \/c r w sm(2n( Zi — Hi)) 
fnA z i) = exp((^ - ni){2ii + h + ej) — , 

2lt(Zi - Hi) 

along with the usual technical provision at each of the removable singular- 
ities. If So is the set of entire functions of C k which vanish at the points 
(zi, . . . , Zk) = (ni, . . . , rik) for < rii G Z and 1 < i < k, the equivalence 
classes have the form S = s + So- Once again they are independent of the 
explicit values of a, bi and used to construct them. The proofs proceed as 
before with only minor alterations and so we omit the explicit details. 



4 Berry's hypothesis 

Comparing the spectral rigidity of quantum systems which are classically 
integrable to those which are classically chaotic, Berry considered the spec- 
trum of an unknown dynamical system whose energy levels are given by the 
imaginary parts of the zeros on the critical line of the Riemann zeta function 
||. It had previously been conjectured by Montgomery that the distribution 
of energy gaps of such a system would be GUE |7| . Montgomery's conjecture 
was supported numerically by the work of Odlyzko according to a report 
by Bohigas and Giannoni ||. Since GUE statistics are normally associated 
with classically chaotic systems which do not possess time reversal invariance, 
this suggested to Berry that the corresponding unknown classical dynamical 
system must be chaotic 0. Berry provided further support for this hypoth- 
esis through theoretical work based on a semi-classical consideration of the 
rigidity of its spectral sequence 0. 
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We will now show that although there may exist classically chaotic sys- 
tems whose quantum spectrum is given by the imaginary parts of the non- 
trivial zeros of the Riemann zeta function, there also exists an infinite family 
of classical integrable systems for which this is true. 

To see that this is so it suffices to show that the monotonic sequence ( n 
where 1/2 + t( n is the n-th nontrivial zero of the Riemann zeta function sat- 
isfies the growth condition of theorem |l[ This assertion follows immediately 
from a classical theorem about the distribution of zeros on the critical line 
due to Hardy and Littlewood ||. Their theorem states that if N(T) is the 
number of zeros of the Riemann zeta function on the interval [1/2, 1/2 + tT], 
then there exists a constant a so that 

N(T) > aT. 

This tells us that if 1/2 + z( n is the position of the n-th zero, then 

C n < a _1 n. 

The sequence {( n :0<n6 Z} is therefore exponentially bounded and satis- 
fies the requirements of theorem [l], which we apply to deduce the following. 

THEOREM 3 There exists an infinite family of classically integrable non- 
linear oscillators whose quantum spectrum is given by the imaginary part of 
the sequence of zeros on the critical line of the Riemann zeta function. 

It is possible to go even further by relabelling the sequence with one 
index ( n as a sequence with k indices Cnx...n k f° r an Y 1 < k E Z , such that 
the growth condition of lemma |2| is still satisfied. There are many ways in 
which to do this and applying theorem |2| we deduce the following. 

THEOREM 4 For any finite dimension k, there exists an infinite family 
of classically integrable k- dimensional non-linear oscillators whose quantum 
spectra reproduce the imaginary part of the zeros on the critical line of the 
Riemann zeta function. 

5 Conclusion 

Apart from an illustration of how the correspondence between classical and 
quantum dynamical systems is many-to-many rather than one-to-many, our 
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main conclusion is that exceptions to the rule of Berry and Tabor regarding 
the distribution of energy gaps in the spectrum of a classically integrable 
system, are more numerous than the literature suggests. In particular we 
show that contrary to the Berry hypothesis, the unknown classical dynamical 
system whose chaotic quantum spectrum is given by the imaginary part of 
the non-trivial zeros of the Riemann zeta function is not unique and need not 
be chaotic. For a given value of Plane's constant and for any finite dimension 
of phase space, there exists an infinite number of classically integrable non- 
linear oscillators whose quantum spectrum simulates that of Berry's unknown 
system. We conclude that the Poisson distribution of energy gaps is not a 
universal property of integrable systems, but of a restricted class of systems 
for which the approximations made by Berry and Tabor are valid. It would be 
of considerable interest to characterise more precisely the range of validity of 
their approximations and consequently of their result regarding the statistical 
distribution of quantum energy gaps of classically integrable systems. 
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